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An angle resolved photoemission study on Ca2Cu02Cl2, a parent compound of high T c supercon- 
ductors is reported. Analysis of the electron occupation probability, n(fc), from the spectra shows a 
steep drop in spectral intensity across a contour that is close to the Fermi surface predicted by the 
band calculation. This analysis reveals a Fermi surface remnant even though Ca2Cu02Cl2 is a Mott 
insulator. The lowest energy peak exhibits a dispersion with approximately the |cos(fc x a)-cos(fc y a)| 
form along this remnant Fermi surface. Together with the data from Dy doped E^S^CaC^Os+a 
these results suggest that this d-wave like dispersion of the insulator is the underlying reason for the 
pseudo gap in the underdoped regime. 



I. INTRODUCTION 

Introduction: A consensus on the d x 2 _ y 2 pairing state 
and the basic phenomenology of the anisotropic normal 
state gap (pseudo gap) in high-T c superconductivity has 
been established |lj , partially on the basis of angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments [|| H Q , in which two energy scales have been 
identified in the pseudo gap, a leading-edge shift of 20-25 
meV and a high-energy hump at 100-200 meV. Both 
of these features have an angular dependence consistent 
with a d-wave gap. For simplicity in the discussion be- 
low, we refer to these as low- and high-energy pseudo 
gaps, respectively, in analogy to the analysis of other 
data. H The evolution of these two pseudo gaps as a 
function of doping are correlated , but the microscopic 
origin of the pseudo gap and its doping dependence are 
still unestablished. Theoretical ideas of the pseudo gap 
range from pre-formed pairs or pair fluctuation j^] and 
damped spin density wave (SDW) || to the evidence of 
the resonating valence bond (RVB) singlet formation and 
spin-charge separation (To) |ll|] p2|| . 

To further differentiate these ideas, it is important to 
understand how the pseudo gap evolves as the doping 
is lowered and the system becomes an insulator. We 
present experimental data from the insulating analog of 
the superconductor La2-xSr x Cu04, Ca2Cu02Cl2 which 
suggest that the high energy pseudo gap is a remnant 
property of the insulator that evolves continuously with 
doping, as first pointed out by Laughlin. The Com- 
pound Ca2Cu02Cl2, a half-filled Mott insulator, has the 
crystal structure of La2CuC>4 and it can be doped by 
replacing Ca with Na or K to become a high-temperature 
superconductor. |l4|] As with the case of S^CuC^C^, 
Ca2Cu02Cl2 has a much better surface property than 
La2Cu04 and thus is better suited for ARPES experi- 
ments. E5[ Although the data from Ca2Cu02Cl2 are con- 
sistent with earlier results from Sr 2 Cu0 2 Cl2, pi |i~7| El 



the improved spectral quality obtainable from this ma- 
terial allows us to establish that: (I) The Fermi surface, 
which is destroyed by the strong Coulomb interactions, 
left a remnant in this insulator with a volume and shape 
similar to what one expects if the strong electron corre- 
lation in this system is turned off; (II) The strong cor- 
relation effect deforms this otherwise iso-energetic con- 
tour (the non-interacting Fermi surface) into a form that 
matches the |cos(fca;a)-cos(A;j,a)| function very well, but 
with a very high energy scale of 320 meV. Thus, a d-wave 
like dispersive behavior exists even in the insulator. 

Comparison with data from underdoped 
Bi 2 Sr2CaCu20 8+(5 (Bi2212) with T c 's of 0, 25 and 65 
K indicates that the high energy d-wave like pseudo gap 
in the underdoped regime originates from the d-wave like 
dispersion in the insulator. Once doped to a metal, the 
chemical potential drops to the maximum of this d-wave 
like function, but the dispersion relation retains its qual- 
itative shape, albeit the magnitude decreases with dop- 
ing. Thus, only the states near the d-wave node touch 
the Fermi level and form small segments of the Fermi sur- 
face, with the rest of Fermi surface gapped. In this way, 
the d-wave high energy pseudo gap in the underdoped 
regime is naturally connected to the properties of the in- 
sulator. Since the high energy pseudo gap correlates with 
the low energy pseudo gap which is likely to be related to 
superconductivity @ [j§ [§ @, it is likely that the 
same physics that controls the d-wave dispersion in the 
insulator is responsible for the d-wave like normal state 
pseudo gap and the superconducting gap in the doped 
superconductors. 



II. METHODOLOGY 

To investigate the strong correlation effect, we contrast 
our experimental data with the conventional results for 
the case when the correlation effects are neglected. Wc 
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can obtain the occupation probability, n(&), by integrat- 
ing K(k,ui) obtained by ARPES, over energy. Q Exper- 
imentally A(k,ui) can not be integrated over all energies 
due to contributions from secondary electrons and other 
electronic states. Instead an energy window for integra- 
tion must be chosen, and the resulting quantity we define 
as the relative n(fc). Fortunately, the features we are in- 
terested in are clearly distinguishable from any other con- 
tributions. We note that n(fc) is a ground state property, 
and hence is different from the integration of the single- 
particle spectral weight, A(fc,w), over energy. However, 
under the sudden approximation integration of A(k,uj) 
as measured by ARPES gives n(fc). (|(]] We then use the 
drop of the relative n(k) to determine the Fermi surface 
as illustrated in Fig. 1. For a metal with non- interacting 
electrons, the electron states are filled up to the Fermi 
momentum, kp, and the n(fc) shows a sudden drop(Fig. 
1A). As more electrons are added, the electron states 
are eventually filled and the system becomes an insula- 
tor with no drop in n(k) (Fig. IB). Therefore, the drop 
in n(k) characterizes the Fermi surface of a metal with 
non-interacting electrons. When correlation increases, 
n(k) begins to deform (Fig. 1C), although there is still 
some discontinuity at k-F when the correlation is moder- 
ate. Note that the electrons that used to occupy states 
below fcj? have moved to the states that were unoccupied. 
For a non-Fermi liquid with very strong correlation, n(k) 
drops smoothly without a discontinuity (panel D). Several 
theoretical calculations using very different models have 
found that n(&) of the interacting system mimics that 
of the non-interacting system, even when the material is 
fully gapped [^l) |2^] Q Q . Hence we can recover the 
remnant of a Fermi surface or an underlying Fermi sur- 
face by following the contour of steepest descent of n(k) 
even when correlation is strong enough that the system 
becomes a Mott insulator. 25 The volume obtained by 
this procedure is consistent with half-filling as expected 
in a Mott insulator. 

We apply this method to determine the Fermi level 
crossing of a real system. The traditional way (Fig. 2A) 
is shown for the ARPES spectra on the (0,0) to (tt, tt) cut 
taken from Bi2212 which is metallic. As we move from 
(0,0) toward (ir, tt), the peak disperses to the Fermi level, 
Fjp. As the peak reaches Ep and passes it, it begins to 
lose spectral weight (this again is Uf)- Alternatively, we 
simply integrate the spectral function from 0.6 eV to -0.1 
eV relative to the Ep, and the resulting relative n(fc) is 
plotted in Fig. 2C. We can now define Uf as the point 
of steepest descent in the relative n(fc). The same con- 
clusion can be drawn here independent of the method we 
use. Note that the n(fc) also drops as we approach (0,0); 
this is a photoemission artifact, because the photoemis- 
sion cross-section of the d x 2_ y 2 orbital vanishes due to 
symmetry. 

We can show that the n(k) procedure is still valid 
for strongly correlated systems with gapped Fermi sur- 



face by presenting ARPES spectra on ferromagnetic 
La^Sr^MnsOT on the (tt, 0) to (7r,7r) cut (Fig. 2B). 
p| It shows a dispersive feature initially moving toward 
Ei? and then pulling slightly back away from Ep around 
(tt, 0.277r), but never reaching the E^. However, the fea- 
ture suddenly loses its spectral weight when it crosses 
(tt, 0.277r) as if it crosses the Fermi surface as shown in 
panel D. Furthermore, the Fermi surface determined by 
a local density approximation calculation coincides with 
the Fermi surface determined by the n(k) despite the 
spectra of this ferromagnetic metallic state material hav- 
ing a significant gap. Thus, the underlying Fermi surface 
can survive a strong interaction, and the n(fc) method is 
effective in identifying it even when the peak does not 
disperse across Ep. 



III. EXPERIMENTAL RESULTS 

The low-energy feature along the (0,0) to (tt, tt) cut 
on Ca2Cu02Cl2 (Fig. 3A) has the same origin as the 
lo- energy peak seen in Bi2212, the Zhang-Rice singlet 
on the CuC>2 plane. As k increases from (0,0) toward 
(tt, tt), the peak moves to lower energy and subsequently 
pulls back to higher energy as it crosses (7r/2,7r/2). Its 
spectral weight increases as it moves away from the (0,0) 
point for the reason described earlier, and then drops as 
it crosses (0.437T, 0.437r). These changes along the (0,0) 
to (tt, tt) cut are consistent with the earlier reports on 
Sr 2 Cu0 2 Cl2. |6| @ Similar to the drop of n(k) 
across the Fermi surface seen in Bi2212, Ca2CuC>2Cl2 
also shows that the intensity of the peak n(k) drops as 
if there is a crossing of Ej? even though the material is 
an insulator. The intensity along the (0,0) to (tt, 0) cut 
(Fig. 3B) goes through a maximum around (27r/3,0) 
as in Sr2Cu02Cl2. This behavior is also seen in super- 
conducting cuprates. Earlier works on Sr2Cu02Cl2 
show the spectral weight along the (tt,0) to (tt,tt) cut is 
strongly suppressed. However, for Ca2Cu02Cl2, the im- 
proved spectral quality allows us to clearly observe the 
spectral weight drop along the (tt,0) to (tt,tt) cut (Fig. 
3C) . [^8) Note that the spectral weight drops as we move 
toward the (tt, tt) point, which we attribute to the cross- 
ing of a remnant Fermi surface. We also show another 
cut (Fig. 3D) which exhibits essentially the same behav- 
ior. The relative n(fc)s of the cuts are summarized in Fig. 
3E in arbitrary units. The relative n(fc) here and in fig. 
4 were obtained by integrating from 0.5 eV to -0.2 eV 
relative to the peak position at (tt/2,tt/2). All of the 
n(fc)s show a drop (after the maximum) as we cross the 
remnant Fermi surface. Here we emphasize that we are 
using the same method as we do for metals, where the 
identification of a Fermi surface is convincing. 

The remnant Fermi surface can be identified in the 
contour plot of n(k) of Ca 2 Cu0 2 Cl 2 (Fig. 4A). The little 
crosses in the figure denote the fc-space points where spec- 
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tra were taken. The data points here and in Fig. 4C have 
been reflected about the line ky=k x to better illustrate 
the remnant Fermi surface. Again, it should be empha- 
sized that the suppressed n(fc) near (0,0) comes from the 
vanishing photoemission cross section due to the d x 2_ y 2 
orbital symmetry rather than a remnant Fermi surface 
crossing. For the same reason, the photon polarization 
suppresses the overall spectral weight along the (0,0) to 
(tt, tt) line as compared with the (0,0) to (tt, 0) line, with 
a monatonic change between the two directions. In fig. 
4B we present the relative n(fc) of an optimally doped 
Bi2212 sample in the normal state. In this case the iden- 
tification of the Fermi surface is unambiguous, but the 
same matrix element effects that were seen in the insula- 
tor can be seen in the metallic sample as well. However, 
for both samples, the drop in n(fc) near the diagonal line 
connecting (tt, 0) and (0, tt) can not be explained by the 
photoemission cross section. In the metallic case, the 
Fermi surface is clearly identified (the white-hashed re- 
gion in Fig. 4B). For the insulator, the drop is approx- 
imately where band theory predicts the Fermi surface. 
p9| Therefore, we attribute the behavior in the insula- 
tor to a remnant of the Fermi surface that existed in the 
metal. The similarity of the results in the insulator and 
the metal makes the identification of the remnant Fermi 
surface unambiguous. The white hashed area in Fig. 4A 
represents the area where the remnant Fermi surface may 
reside as determined by the relative n(k). Although there 
is some uncertainty in the detailed shape of this remnant 
Fermi surface, this does not affect the discussion and the 
conclusions drawn below. The relative n(fc) we presented 
is a very robust feature. In metallic samples with par- 
tially gapped Fermi surfaces, underlying Fermi surfaces 
have also been identified in the gapped region. |^6| ]3C| ] 
This effect is similar to what we report here in the in- 
sulator. The remnant Fermi surface in the underdoped 
Bi22I2 was also identified at similar locations to the n(k) 
drop in these materials with a different criteria of min- 
imum gap locus. (3^] Calculations also show the Fermi 
surface defined by n(k) is robust in the presence of strong 
correlation. j2l| |23) Given that there is a remnant 
Fermi Surface as shown by the white hashed lines in Fig. 
4, A and C, the observed energy dispersion along this 
line has to stem from the strong electron correlation. In 
other words, the electron correlation disperses the other- 
wise iso-energetic contour of the remnant Fermi surface. 
This dispersion is consistent with the non-trivial rf-wave 
|cos(A: :I ;a)-cos(/; y a)| form. 10 |l2| These results also 
support our identification of the remnant Fermi surface 
in a Mott insulator. 

Fig. 4C plots the energy contour of the peak position 
of the lowest energy feature of Ca2CuC>2Cl2 referenced 
to the energy of (tt/2, tt/2) peak. The hashed area indi- 
cates the remnant Fermi surface determined in fig. 4A. 
The 'Fermi surface' is no longer a constant energy con- 
tour as it would be in the non-interacting case. Instead 



it disperses as much as the total dispersion width of the 
system. In Fig. 4D we plot the dispersion at differ- 
ent points on the remnant Fermi surface referenced to 
the lowest energy state at (tt/2, tt/2). The dispersion 
of the peaks along the Fermi surface is plotted against 
cos(A; a ;a)-cos(A:j,a)|. The straight line shows the d-wave 
dispersion function at the 'Fermi surface' with a d-wave 
energy gap. The figure in the inset presents the same 
data in a more illustrative fashion. On a line drawn from 
the center of the Brillouin zone to any point either ex- 
perimental (blue) or theoretical (red), the distance from 
this point to the intersection of the line with the anti- 
ferromagnetic Brillouin zone boundary gives the value of 
the 'gap' at the fc-point of interest. The red line is for a 
d-wave dispersion along (7r,0) to (0,7r). The good agree- 
ment |3l| l is achieved without the need for free parame- 
ters. This d-wave like dispersion can only be attributed 
to the many-body effect. The relative energy difference 
between the energy at (7r/2, tt/2) and (tt, 0) has been re- 
ferred to as a gap [Q, which we follow. 

This gap differs from the usual optical Mott gap (Fig. 
5) and may correspond to the momentum dependent gap 
once the system is doped. This gap monotonically in- 
creases when we move away from (tt/2, tt/2) as also re- 
ported earlier. |16| As well as summarizing the data pre- 
sented, Fig. 5 also shows the intriguing similarity be- 
tween the data from the insulator and a slightly over- 
doped d-wave superconductor(Bi22I2), and thus gives 
the reason for comparing the dispersion along the rem- 
nant Fermi surface with the |cos(fc :r a)-cos(fc y a)| form. In 
the superconducting case, n(k) helps determine the Fermi 
surface. The anisotropic gapping of this surface below T c 
reveals the d-wave nature of the gap. In the insulator, 
n(fc) helps determine the remnant Fermi surface. The 
fc-dependent modulation along this surface reveals the d- 
wave like dispersion. Whether this similarity between the 
insulator and the doped superconductor is a reflection of 
some underlying symmetry principle is a question which 
needs to be investigated. [g2| 

The above analysis is possible only because we now 
observe the remnant Fermi surface. Although the dis- 
persion for Sr2Cu02Cl2 was similar to the present case, 
the earlier results did not address the issue of a remnant 
Fermi surface because the smaller photoemission cross 
section along the (tt, 0) to (tt, tt) cut prevented this iden- 
tification. Therefore the analysis shown above was not 
possible. With only the energy contour information (such 
as in Fig 4C), it is plausible to think that the Fermi sur- 
face evolves to a small circle around the (tt/2, tt/2) point. 
p3| However, with the favorable photoemission cross sec- 
tion, the results from Ca2Cu02Cl2 show that the Fermi 
surface leaves a clear remnant, although it may be broad- 
ened and weakened. Therefore, the energy dispersion 
along the original Fermi surface of a non-interacting sys- 
tem is due to the opening of an anisotropic 'gap' along 
the same remnant Fermi surface. 
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The same analysis is shown in Fig. 6 for Bi2212 
with different Dy dopings together with Ca2Cu02Cl2 re- 
sults. The corresponding doping level and T c as a func- 
tion of Dy concentration are also shown. The energy 
for Ca2Cu02Ci2 is referenced to the peak position at 
(tt/2, tt/2) and that for Dy doped Bi2212 is to E F . How- 
ever, the two energies essentially refer to the same energy 
since the peak on the (0,0) to (tt,tt) cut for all Bi2212 
samples reaches the Fermi level. Note that the gaps for 
Dy doped Bi2212 data also follow a function that is qual- 
itatively similar to the d-wave function with reduced gap 
sizes as shown with the (tt, 0) spectra in Fig. 6B. This 
result suggests that the d-wave gap originating in the in- 
sulator continuously evolves with doping, but retains its 
anisotropy as a function of momentum and that the high 
energy pseudo gap in the underdoped regime is the same 
gap as the d-wave gap seen in the insulator as discussed 
above. Of course, the high energy pseudo gap in the un- 
derdoped regime is smaller than the gap in the insulator. 
In a sense, the doped regime is a diluted version of the 
insulator, with the gap getting smaller with increasing 
doping. The two extremes of this evolution are illus- 
trated in the quasiparticle dispersions shown in Fig. 5. 
The insulator shows a large d-wave like dispersion along 
the remnant Fermi surface. In the overdoped case, no 
gap is seen in the normal state along an almost identi- 
cal curve in A;-space; however, a d-wave gap is observed in 
the superconducting state. Although their sizes vary, the 
d-wave superconducting gap, and the d-wave 'gap' of the 
insulator have the same non-trivial form, and are thus 
likely to stem from the same underlying mechanism. 

IV. DISCUSSION 

We do not know the full implications of the data we 
report, but can offer the following possibilities. First, we 
compare the experimental dispersion with a simple spin- 
density wave picture. Starting with the Hubbard model 

H = Sfco- €k c' kaCkcr + US^ n^n^ 

with 

efc = -2t(cos(fca;a)-l-cos(A;j,a)) - 4t'(cos(fc a; a)cos(fcya)) - 
2t" (cos(2k IE a)+cos(2kj / a)) 

and adding a SDW picture, the following dispersion 
relation will be found 

E fc± 

« -4t'(cos(fc2;a)cos(fcj,a)) - 2t" (cos(2k a; a)-|-cos(2kj / a)) ± 
[U/2 + 3(cos(k x a,)+cos(k y a,)) 2 } 

with J=t 2 /U. With realistic values for t' and t", and 
an experimental value for J of -0.12 eV, 0.08 eV, and 
0.125 eV respectively, Jl8| we find that the experimental 
dispersion deviates significantly from this mean field re- 
sult giving a bandwidth of 1.1 eV. It is crucial to note 
the observed isotropic dispersion around the (tt/2, tt/2) 
point, with almost identical dispersions from (tt/2, tt/2) 



to (0,0) and from (7r/2,7r/2) to (tt, 0). This result is un- 
likely to be a coincidence of the parameters t', t" , and J 
as suggested by the SDW picture above. 

We now compare the data with numerical calculations 
that, unlike the mean field SDW picture, appropriately 
accounts for the dynamics. Being mainly concerned with 
the dispersion relation, we concentrate our discussion on 
the t-J model as more extensive literature exists and as 
J can be independently measured. |^4| Qualitatively, the 
same conclusion is expected for the Hubbard model |34| , 
which has the added advantage of yielding n(k), but has 
more uncertainty in the parameter U. Although the t-J 
model correctly predicts J24j the dispersion along (0,0) 
to (7r/2,7r/2) quantitatively [pUf , with the band width 
along this direction solely determined by J, it incorrectly 
predicts the energy of (tt, 0) to be nearly degenerate to 
(tt/2, tt/2). This is a serious deficiency of the t-J model, 
because the evolution of the (tt,Q) feature is crucial to 
understand the d-wave-like pseudo gap. The inclusion 
of the next nearest neighbor hoppings of t' and t" can 
resolve this problem. @ || |(§ In fact, the t-t'-t"-J 
model can account for both the dispersion and lineshape 
evolution over all doping levels, which is a remarkable 
success of this model. |l^] |3j| With a J/t ratio in the 
realistic range of 0.2 to 0.6, the t-t'-t"-J model shows 
that the dispersion from (tt/2, tt/2) to (0,0) and to (tt,0) 
are equal and scaled by J. |37j This result supports the 
notion that the isotropic dispersion is controlled by a 
single parameter, J, as stressed by Laughlin. Jl2| 

The above discussion indicates that we have a model, 
when solved by Monte Carlo or exact diagonalization, 
that can account for the data, but what does the data 
fitting the non-trivial |cos(fc 3; a)-cos(fcj / a)| function so well 
mean? As pointed out [Q, the key to the inclusion of 
t' and t" is that the additional hole mobility destabilizes 
the one-hole Neel state with the hole at (tt, 0) and makes 
the system with one-hole move closer to a spin liquid 
state rather than to a Neel state that is stable in the t-J 
model. This point is relevant to some early literature of 
the resonating valence bond(RVB) state j3f| |3j|. Ander- 
son conjectured that the ground state of the insulator at 
half filling is a RVB spin liquid state. |3^] This idea was 
extended in the context of a mean field approach to the t- 
J model that yields a d-wave RVB or flux phase solution. 
ifiof The mean-field solution also predicts a phase dia- 
gram similar to what is now known about the cuprates, 
with the d-wave like spin gap in the underdoped regime 
being the most successful example. The problem with 
the mean-field solution of the t-J model is that it docs 
not agree with exact numerical calculation results pi] , 
and the half-filled state was found by neutron scattering 
to have long range order. |^] If these numerical calcu- 
lations are right then the d-wave RVB is not the right 
solution of the t-J model. However, the d-wave RVB like 
state may still be a reasonable way to think about the 
experimental data that describes the situation of the spin 
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state near a hole. jy] It is just that one has to start with 
a model where the single hole Neel state is destabilized, 
as in the t-t'-t"-J model. We leave this open question as 
a challenge to theory. 

The presence of c?-wave like dispersion along the rem- 
nant Fermi surface shows that the high energy pseudo 
gap is a remnant of the rf-wave 'gap' seen in the insu- 
lator. The details of the evolution of this gap, and its 
connection to the low energy pseudo gap (which is likely 
due to pairing fluctuations) as well as the superconduct- 
ing gap is unclear at the moment. However, we believe 
that there has to be a connection between these gaps of 
the similar |cos(fca;a)-cos(A;j,aj| form, as their presence is 
correlated with each other. 
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FIG. 1. Illustration of the Fermi surface determination. 

(A) The case for band metal. Electrons occupy states only 
up to a certain momentum, showing a sharp drop in n(fc). 

(B) Band insulator case. Electrons occupy all possible states 
and do not show a drop in n(k). (C) Fermi liquid with elec- 
tron correlation. Note that electrons that used to occupy the 
states below fcj- have moved above fcf- However, it still shows 
a discontinuity at &f- (D) For a strongly correlated non- Fermi 
liquid n(fc) does not show discontinuity, yet there exists n(k) 
drop showing the remnant Fermi surface. 
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FIG. 2. Application of the method described in Fig.l. (A) 
Spectra along the (0,0) to (7r,7r) cut from Bi2212. The peak 
disperses towards the low energy side and reaches the Fermi 
level at k,F, (0.437T, 0.437r). (B) Spectra along the (tt,0) to 
(n, 7r) cut from metallic Las-^Sr^M^Oy. The peak disperses 
toward the low energy side, but never reaches the Fermi en- 
ergy. However, it loses intensity as it crosses the position 
where the band calculation predicts the Fermi surface, show- 
ing an underlying Fermi surface. (C) and (D) plots of the 
relative n(k) for the data in (A) and (B), respectively, show a 
sudden drop around Uf , essentially showing the two methods 
give the same Fermi momentum Icf- 




FIG. 3. ARPES Spectra on various cuts from Ca2Cu02Ci2 
and n(k) plots. The insets and labels show where the spectra 
were taken in the Brillouin zone. (A) (0,0) to (it, it) cut. The 
peak disperses towards the low energy side and loses intensity 
near the (tt/2,it/2) point. (B) (0,0) to (it,0) cut. The lowest 
energy peak shows little dispersion. The spectral weight ini- 
tially increases and then decreases again after (0.677T, 0) as in 
the Sr2Cu02Cl2 case. [17] [18] However, note that there is ap- 
preciable spectral weight at (it, 0) contrary to the Sr2Cu02Ci2 
case. (C) (tt,0) to (it, it) cut. The spectral weight drops as 
we move to (it, it). (D) Another cut (as marked in the inset) 
showing the n(k) drop. (E) relative n(fc)'s constructed from 
the spectra in panels A-D. The relative increase of spectral 
weight above is caused by emission from second order light. 
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FIG. 4. Contour plot of the relative n(fc). (A) n(k) from the 
spectra shown in fig. 3. The color scale on the right represent 
n(k). The spectra were taken only in the first octant of the 
first Brillouin zone (crosses). The n(k) plot was folded to 
better represent the remnant Fermi surface. Note, the n(k) 
drops as we cross the approximate diagonal line connecting 
(0, 7r) and (n,0). The hashed area represents approximately 
where the remnant Fermi surface exists. (B) An identical plot 
for an optimally doped Bi22I2 sample in the normal state. 
(C) Contour plot of the lowest energy peak position from the 
spectra in fig. 3 relative to the (n/2,n/2) peak. The hashed 
area is from (A), showing the remnant Fermi surface. The 
color scale on the right indicates the relative binding energy 
of the peak. The peak disperses isotropically away from the 
(tt/2, 7r/2) peak position as we move away from the (7r/2, 7t/2) 
point. (D) The 'gap' versus |cos(fc a; a)-cos(i;j / a)|. The straight 
line shows the d-wave line. The inset is a more illustrative 
figure of the same data as explained in the text. 



FIG. 5. An illustration showing the 2 experimental features 
presented in this paper on the insulator, and the similarity 
they show to a slightly overdoped Bi22I2 sample. (A) The 
bottom half shows the relative n(A;), and above it lies the 
approximate remnant Fermi surface derived from it. However, 
there is much dispersion over the entire Brillouin zone, and the 
remnant Fermi surface is no longer an iso-energetic contour as 
can be seen by the quasiparticle dispersion (energy relative to 
the valence band maximum) . Here the remnant Fermi surface 
is shown as a black and white line running over the visible 
portions of the dispersion contour. For clarity, a portion of 
the dispersion along the remnant Fermi surface is shown in 
the top half. Note the idea presented that the isoenergetic 
contour(dashed black line) is deformed by strong correlation 
to the observed red curve. The d-wave like 'gap' referred 
to in the text is the quasiparticle energy deviation from the 
dashed black line set at the energy of the (7r/2,7r/2) point. 
The difference between this 'gap' and the Mott gap can now 
be seen clearly. (B) For overdoped Bi22I2, n(fc) defines the 
actual Fermi surface. The quasiparticle dispersion (binding 
energy) shows states filled to an isoenergetic Fermi surface. 
In the top panel one is reminded that below T c , a d-wave 
superconducting gap opens. This is an intriguing similarity 
between the insulator and the metal. 
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FIG. 6. (A) Combined a!-wave plot of the data from 
Ca 2 Cu0 2 Cl2 and Bi2212 with various Dy dopings. (B) The 
spectra at (-7T, 0), showing the evolution of the high energy 
pseudo gap as a function of doping, as previously stressed by 
Laughlin. [12] 
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